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Abstract. It is by now well known that the wave functions of rational solu- 
tions to the KP hierarchy which can be achieved as limits of the pure n-soliton 
solutions satisfy an eigenvalue equation for ordinary differential operators in 
, the spectral parameter. This property is known as "bispectrality" and has 

Q i proved to be both interesting and useful. In this note, it is shown that all pure 

►7- ' soliton solutions of the KP hierarchy (as well as their rational degenerations) 

satisfy an eigenvalue equation for a non-local operator constructed by compos- 
ing ordinary differential operators in the spectral parameter with translation 
operators in the spectral parameter, and therefore have a form of bispectrality 
as well. 

> 

On 
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1. Introduction 

1.1. The KP Hierarchy and Bispectrality. Let D be the vector space spanned 
over C by the set 

{A(j,A) I AeC,jeN} 

whose elements differentiate and evaluate functions of the variable z: 

A(j, A)[/(z)] :=/«(A). 

O ' The elements of B are thus finitely supported distributions on appropriate spaces of 

functions in z. For lack of a better term, we will continue to call them distributions 
even though their main use in this paper will be their application to functions of 
two variables. (Such distributions were called "conditions" in |20| since a KP wave 
\— 1 ' function was specified by requiring that it be in their kernel.) Note that if c s B and 

fix, z) is sufficiently differentiable in z on the support of c, then fix) — c[f(x, z)] 
is a function of x alone. Furthermore, note that one may "compose" a distribution 
with a function of z, i.e. given c G B and f(z) (sufficiently differentiable on the 
support of c) then there exists a d := c o / s B such that 

c'{g(z)) = c(fiz)g(z)) V.g. 

The subspaces of B can be used to generate solutions to the KP hierarchy in 
the following way. Let C C B be an n dimensional subspace with basis {or, ... , c„}. 
Then, if K = Kq is the unique, monic ordinary differential operator in x of order 
n having the functions Ci(e xz ) in its kernel (see (|J)) we define Cc = K^K -1 
and ipc = -^Ke xz . The connection to integrable systems comes from the fact that 
adding dependence to C on a sequence of variables tj (J = 1,2,...) by letting C(tj) 
be the space with basis 

{ci o e S ~ t]Zl , c 2 o e S ~ tjzJ , c„ o e S ~ t]Z ° } 
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it follows that the "time dependent" pseudo-differential operator C = C(tj) satisfies 
the equations of the KP hierarchy [ll], O, |2(| 

9 -£ = {(&)+,£}■ 



atj 



The wave function tpc{x, z) generates the corresponding subspace of the infinite di- 
mensional grassmannian Gr |L8| which parametrizes KP solutions and thus it is not 
difficult to see that this construction produces precisely those solutions associated 



to the subgrassmannian Gr\ C Gr 18, [20| . 

Moreover, the ring Ac = {p E C[z]\aop s C 1 < i < n} is necessarily non-trivial 
(i.e. contains non-constant polynomials) and the operator L p = p(£) is an ordinary 
differential operator for every p G Ac and satisfies 

(1) Lptpcix, z) = p(z)ipc(x, z). 

The subject of this paper is the existence of additional eigenvalue equations satisfied 
by il>c(x, z). In particular, we wish to consider the question of whether there exists 
an operator A acting on functions of the variable z such that 

(2) Atpc{x,z) = n(x)tpc(x,z) 

where tt(x) is a non-constant function of x. For example, the following theorem is 
due to G. Wilson in pp|: 

Theorem 1.1. In addition to ([!]) the wave function ipc(x, z) is also an eigenfunc- 
tion for a ring of ordinary differential operators in z with eigenvalues depending 
polynomially on x if and only if C has a basis of distributions each of which is 
supported only at one point. 

In other words, for this special class of KP solutions for which the coefficients of £ 
are rational functions of x, the wave function ipc satisfies an additional eigenvalue 
equation of the form (Q) where A is an ordinary differential operator in z and 
7r(x) a non-constant polynomial in x^ Together ([!]) and (0) are an example of 
bispectrality fi |j. The bispectral property is already known to be connected to 
other questions of physical significance such as the time-band limiting problem in 
tomography [pj, Huygens' principle of wave propagation [0], quantum integrability 
|)| [T^] and, especially in the case described above, the self duality of the Calogero- 
Moser particle system |(], |l| . 

It is known that the only subspaces C for which the corresponding wave func- 
tion satisfies (|l|) and (|^) with L p and A ordinary differential operators in x and z 



respectively are those described in Theorem 1.1. However, suppose we allow A to 
involve not only differentiation and multiplication in z but also translation in z and 
call this more general situation t-bispectrality.^j It will be shown below that there 
are more KP solutions which are bispectral in this sense. In particular, it will be 
shown that the KP solution associated to any subspace C shares its eigenfunction 
with a ring of translational-differential operators in the spectral parameter. 



1 Moreover, he demonstrated that up to conjugation or change of variables, the operators L p 
found in this way are the only bispectral operators which commute with differential operators of 
relatively prime order, but this fact will not play an important role in the present note. 

2 It should be noted that the term "bispectrality" already applies to more general situations 
than simply differential operators Q , but in the case of the KP hierarchy I believe only differential 
bispectrality has thus far been considered. 
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1.2. Notation. Using the shorthand notation d = ^ any ordinary differential 
operator in x can be written as 

N 

L = Y,f^W (NeN). 

i=0 

We say that a function of the form 

n 

f(x) = Y,Pi(x)e XiX Aj G C, ^ G C[ar] 
i=i 

is a polynomial- exponential function and that the quotient of two such functions 
is rational- exponential. This note will be especially concerned with the ring of 
differential operators with rational-exponential coefficients and especially with the 
subring having polynomial-exponential coefficients. Similarly, we will write d z = 
but will need to consider only differential operators in z with rational coefficients. 

For any A G C let S> = e xdz be the translational operator acting on functions of 
z as 

S x [f(z)} = f(z + X). 

Then consider the ring of translational-differential operators T generated by these 
translational operators and ordinary differential operators in z. Any translational- 
differential operator T G T can be written as 

JV 

f = 5>i(z,0*)S Ai 

where pi are ordinary differential operators in z with rational coefficients and 
N G N. Note that the ring of ordinary differential operators in z with rational 
coefficients is simply the subring of T of all elements which can be written as pS 
for a differential operator p. 

2. Translational Bispectrality of C[d] 

It has been frequently observed that the ring A = C[d] of constant coefficient 
differential operators in x is bispectral since it has the eigenfunction e xz which it 
shares with the ring of constant coefficient differential operators in z. Here, however, 
we will consider a more general form of bispectrality for the ring A. 

Let A' C T be the ring of constant coefficicient translational-di&crcntisl opera- 
tors. Note that for any element T G A' of the form 

N 
i=l 

one has simply that 

f[e xz ] = (^2 Pt (x)e x * x ^J e xz . 

In particular, e xz is an eigenfunction for the operator with an eigenvalue which is a 
polynomial-exponential function of x. Consequently, the rings A and A 1 are both 
bispectral, sharing the common eigenfunction e xz . 

Let 1Z be the ring of differential operators in x with polynomial-exponential co- 
efficients and TV be the ring of translational-differential operators in z with rational 
coefficients. Note that 1Z is generated by A and the eigenvalues of the operators in 
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A 1 while TZ' is generated by A' and the eigenvalues of the elements of A. It then 
follows @ (see also |]§) that the map b : TZ — > TZ' defined by the relationship 

L[e xz ] = b{L)[e xz ] VL £ K 

is an anti- isomorphism of the two rings. 

3. Translational Bispectrality of KP Solitons 

Let us say that a finite dimensional subspace C C D is t-bispectral if there 
exists a translational-differential operator A 6 T satisfying equation (0) for the 



corresponding KP wave function ipc{x,z). By Theorem 1.1 and the fact that the 
ring of rational coefficient ordinary differential operators in z is contained in T, we 
know that C is t-bispectralQ if it has a basis of point supported distributions. Here 
we will show that, in fact, all subspaces Ccl are f-bispectral. 

An important object in much of the literature on integrable systems is the "tau 
function". The tau function of the KP solution associated to C can be written 
easily in terms of the basis {c{\. In particular, define (cf. |20| ) 

T C (x)=Wr( Cl (e xz ),c 2 (e xz ),... .^(e**)) 

to be the Wronskian determinant of the functions Ci(e xz ). Similarly, there is a 
Wronskian formula for the coefficients of the operator Kq since its action on an 
arbitrary function f(x) is given as: 

(3) Kc(f(x)) = -^m{c 1 (e?%c 2 {<?%...,c n (e?'),f(x)). 

T C (X) 

Then the coefficients of the differential operator Kc ■— Tc(x)Kc(x,d) are all 
polynomials-exponential functions . 

Lemma 3.1. For any C C D there exists a constant coefficient operator L £ A 
which factors as 

1 

where Q g , Kc £ TZ and n(x) = g(x)rc{x) 6 TZ is a polynomial- exponential function. 



L =QgO —f-f o K c 



Proof. Let Aj G C (1 < i < N) be the support of the distributions in C and mi be 
the highest derivative taken at A^ by any element of C. Then the polynomial 

(4) qc(z) := (z - A 4 ) m ' +1 

has the property that c o qc = for any c £ C. Let L$ :— qc(d) and consider 
Lo[c(e xz )] for any element c £ C. Since Lq is an operator in x alone, it commutes 
with c and we have 

L [c(e x *)} =c(L [e xz ]) = c(q(z)e**) =coq(e xz ) = 0. 

So, by the definition of Kq, we see that Lq annihilates the kernel of Kc and thus has 
a right factor of Kc- This gives a factorization of the form Lq = QoKc with Q hav- 
ing rational-exponential coefficients. Then, by choosing a polynomial-exponential 
function g(x) so that Q g := Q o g(x) £ TZ we find the desired factorization. □ 

Given this factorization, the t-bispectrality of all C's now follows from Theorem 
4.2 in @. 

3 ...and also bispectral in the sense of |20|. 
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Theorem 3.1. For any subspace C C D the corresponding KP wave function 
%pc{x,z) satisfies the eigenvalue equation 

^g[i>c{x,z)] = g(x)r c (x)ip c (x,z) 

where A g G T is the translational- differential operator defined by 

n 

A g := z- n o b(K c ) o b(Q g ) o - 



qc{z) 



with Q g defined as in Lemma 3.1. 



Proof. Formally introducing inverses we determine from Lemma 3.1 that 

n(. x ) 9{x)tc{x) = K c o Lq 1 o Q 
and hence (by applying the anti-involution b to this equation) 

b(7r(x)) =b(Q)o-L-ob(K c ). 
Qc(z) 



Then 



z n 1 

k a tyc{x,z)\ = z- n ob(Kc)ob(Q)o—-[^—-Kce x *} 

qc(z) z n Tc(x) 

z~ n - 1 

■ ob(K c )ob(Q)o— -[K c e xz ] 



Tc{x) qc{z) 
ob(Kc)[n(x)e x 



T C (x) 

z- n ir(x) 



tc(x) 

7T(x)tjjc(x, z) 



K c [e 



□ 



Note that according to Theorem [3.1| , each operator A g satisfies an intertwining 
relationship 

Wob(ir(x)) = AgoW 

with the constant coefficient operator 6(7r(a;)) where W = z~ n o b(Kc). As a result 
we find that: 

Corollary 3.1. The set of all such operators A g for a given subspace CcB form 
a commutative ring of translational- differential operators. 

4. Examples 

If we choose C to be the two dimensional space spanned by c\ — A(1,0) and 
C2 = A(I, I) (a "two-particle" Calogero-Moser type solution) then 

, . , 2 + x - (2x + x 2 )z , TZ 
i/>c(x, z) = (l + i-= >-)e xz . 



In this case the translational-differential operators A given by Theorem 3.1 are 
simply ordinary differential operator. For instance, 

l Q q 3 n2 6z 2 -12z + 3 n 122-6 
A = di H -dt —, -T^d 



Z3(Z-1) 2 Z 2 { Z -lf 
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which satisfies Kipc(x, z) = x ipc(x, z) (as we would expect from earlier results on 
bispectrality.) 

However, if we had chosen instead c\ = A(0, 1) + A(0, —1) and C2 = A(0, 2) + 
A(0, 0) we would instead have the case of a 2-soliton solution with 

, / x h 6 + (3z - 2)e 2x + 2z - ze- 2x . xz 

^{x,z) = (i (e , +e -, )2z2 — > xz - 

One finds from the procedure given in the theorem that 

A = z- 2 o ((20z + llz 2 -8z 3 + z 4 )S_ 3 + (60 - 68z - z 2 + 8z 3 + z 4 ) S 5 

+ (-36 + 24z + 16z 2 - 16z 3 + 4z 4 ) S_i + (-44 - 88z - 8z 2 + 16z 3 + 4z 4 ) S 3 

(-12 - 16, -2, 2 + 6* 4 )s0 V _ 2 /!, 2 + 2 , 
satisfies A[ipc(%, z)] = e~ 3x (l + e 2a: ) 4 V>c(£, z). 



5. Conclusions 

In addition to being a generalization of the results of || ^| on bispectral ordinary 
differential operators, the present note may be seen as a generalization of fl5f| in 
which wave functions of n-soliton solutions of the KdV equation are shown to satisfy 
difference equations in the spectral parameter. The idea that KP solitons might be 
translationally bispectral was proposed in [jl4| . 

As in H |2Cj, the equations (Q) and (||) lead to the well known "ad" relations 
associated to bispectral pairs. That is, defining the ordinary differential operator 
A m in x and the translational-differential operator A m in z by 

A m = adZ>(*)) A m = (-lrad^A) 

one finds that A m ij}c{x, z) = A m ipc(x, z). Similarly, if 

B m = ad™ (x) (i p ) B m = (-l) m ad™(p(z)) 

then B m t/jc(x, z) = B m tjjc(x, z). Note that whenever the order of S m _i = N > 
the order of B m cannot be greater than N — 1. So, the familiar result that B m = 
and B m = for m > ordL p holds, which is clearly a strong restriction on the 
operator A. However, unlike the case of bispectral ordinary differential operators, 
one cannot conclude that A m = for sufficiently large m since the order of A m 
may not be reduced by increasing m. 

The bispectrality of the rational KP solutions |20) has been shown to have a 
dynamical significance. In particular, it was shown that the bispectral involution 
is the linearizing map for the classical Calogero-Moser particle system jy], ^lj . 
Moreover, other bispectral KP solutions have been found to have similar proper- 
ties H 0. This would seem to indicate that it is likely that the bispectrality of 
KP solitons also has a dynamical significance, as a map between the classical Rui- 
jsenaars and Sutherland systems (cf. p7| ). In fact, such a bispectral relationship 
between the quantum versions of these systems has been recently found in || . The 
dynamical significance of these results will be considered in a separate paper. 



SPECTRAL DIFFERENCE EQUATIONS SATISFIED BY KP SOLITON WAVEFUNCTIONS 7 



References 

B. Bakalov, E. Horozov and M. Yakimov, "General methods for constructing bispectral op- 
erators", Phys. Lett. A 222 (1996), no. 1-2, 59-66. 

B. Bakalov, E. Horozov and M. Yakimov, "Backlund-Darboux Transformations in Sato's 
Grassmannian" , Serdica Math. Journal, 22 (1996) pp. 571-588 

B. Bakalov, E. Horozov and M. Yakimov, "Commutative Rings of Bispectral Ordinary Dif- 
ferential Operators", to appear Comm. Math. Phys, (1998). 

Yu. Berest, "Huygens' principle and the bispectral problem" , CRM Proceedings and Lecture 
Notes, volume 14, American Mathematical Society, Providence, RI, pp. 11-30 
J.J. Duistermaat and F.A. Griinbaum, "Differential equations in the spectral parameter", 
Comm. Math. Phys. 103 (1986) 177-240 

O.A. Chalykh, "Duality of the generalized Calogero and Ruijsenaars problems" Uspekhi Mat. 
Nauk 52 (1997), no. 6(318), 191-192 

F.A. Griinbaum, "Time-band limiting and the bispectral problem" , Comm. Pure Appl. Math 
47 (1994) no3. 307-328. 

F. A. Griinbaum, "Bispectral Musings" CRM Proceedings and Lecture Notes, volume 14, 
American Mathematical Society, Providence, RI, 1988 pp. 31-46. 

E. Horozov and A. Kasman, "Darboux Transformations for Bispectral Quantum Intcgrablc 
Systems" , in preparation. 

A. Kasman, "Rank r KP Solutions with Singular Rational Spectral Curves", Ph.D. thesis, 
Boston University (1995) 

A. Kasman, "Bispectral KP Solutions and Linearization of Calogero-Moser Particle Systems" , 
Communications in Mathematical Physics 172 (1995) pp. 427-448 

A. Kasman, "The Bispectral Involution as a Linearizing Map", to appear, Proceedings of 
Workshop on Calogero-Moser-Sutherland Models, J.F. van Diejen and L. Vinet eds., Springer- 
Verlag (1998). 

A. Kasman and M. Rothstcin, "Bispectral darboux transformations: the generalized Airy 
case", Physica D 102 (1997) p. 159-173. 

A. Kasman, "Bispcctrality and KP Solitons", Preprint CRM-2533 

M. Reach "Difference equations for AT-soliton solutions to KdV" Phys. Lett. A 129 (1988), 
no. 2, pp. 101-105 

M. Rothstein, "Explicit Formulas for the Airy and Bessel Bispectral Involutions in Terms of 
Calogero-Moser Pairs", CRM Proc. and Lecture Notes, Vol. 14, (1998) pp. 105-110 
S.N.M. Ruijsenaars, "Action-Angle Maps and Scaterring Theory for Some Finite-Dimensional 
Integrable Systems", Commun. Math. Phys. 115 (1988) pp. 127-165 

G. Segal and G. Wilson, "Loop groups and equations of KdV type" Publications Mathema- 
tiques No. 61 de I'Institut des Hautes Etudes Scientifiques (1985) pp. 5-65 

A. P. Veselov, Baker-Akhiczer functions and the bispectral problem in many dimensions, CRM 
Proceedings and Lecture Notes, volume 14, American Mathematical Society, Providence, RI, 
pp. 123-129. 

G. Wilson, Bispectral commutative ordinary differential operators, J. Reine Angew. Math. 
442 (1993), 177-204. 

G. Wilson, Collisions of Calogero-Moser Particles and an Adelic Grassmannian, Inventiones 
Mathematicae 133 (1998), no. 1, 1-41 



Mathematical Sciences Research Institute, Berkeley, CA 94720 



